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Superconductivity induced structures in the electronic Raman spectra of high-T c superconductors 
are computed using the results of ab initio LDA-LMTO three-dimensional band structure calcula- 
tions via numerical integrations of the mass fluctuations, either in the whole 3D Brillouin zone or 
limiting the integrations to the Fermi surface. The results of both calculations are rather simi- 
lar, the Brillouin zone integration yielding additional weak structures related to the extended van 
Hove singularities. Similar calculations have been performed for the normal state of these high-T c 
cuprates. Polarization configurations have been investigated and the results have been compared to 
experimental spectra. The assumption of a simple d x 2_ y 2-like gap function allows us to explain a 
number of experimental features but is hard to reconcile with the relative positions of the Ai g and 
Big peaks. 



I. INTRODUCTION 



After ten yeaxs of research in the field of high-T c 
superconductorsu (HTSC), many of their properties have 
not yet been understood. rip. particular, the symmetry of 
the superconducting gapcm is still controversial. Usu- 
ally, one assumes that the superconducting condensate 
in the HTSC can be described by an order parameter 
Ak, which depends only on the quasi momentum k, but 
not on band index n. Retardation effects are also often 
neglected, i.e. the gap is assumed to be independent of 
the frequency u>. 

A wide range of experimental techniques can be em- 
ployed to investigate the properties of the gap function. 
Among these, Raman scattering has played an important 
role.crQ The dependence of the Raman response on the 
directions of polarization of the incident and scattered 
light yields several independent spectra which provide a 
considerable number of constraints on the assumed k- 
dependence of the gap function Ak- However, Raman 
scattering is not sensitive to the phase of the gap. 

The Raman spectra at temperatures below T c shows, 
in most HTSC, a clear gap-like structure which lies in 
the energy range of the optical phonons at the r point. 
These phonons have been identified for most HTSC,Q and 
the subtraction of the corresponding structures from the 
spectra has become a standard procedure to isolate elec- 
tronic structures containing gap information. Electronic 
Raman scattering spectra are now available for many 
high-T c materials and, since they exhibit similar general 
features, most of these data are considered to be reli- 
able. In this paper, we attempt to interpret these spec- 
tra from a theoretical point of view based on the full 3D 
one-electron band structure. We pay attention to both, 
line shapes and absolute scattering efficiencies. 

The theory of electronic Raman scattering in super- 



conductors was pioneered by Abrikosov and coworkers in 
two important papers flu In the first, they developed a 
theory for the scattering efficiency of isotropic Fermi liq- 
uids under the assumption that the attractive interaction 
between quasiparticles can be neglected. In the second 
paper, they extended this approach to anisotropic sys- 
tems, introduced the effective mass vertex concept, and 
included Coulomb screening. The current form of the 
theory, developed mainly by Klein et aL,£2l takes into ac- 
count the attractive pairing interaction and emphasizes 
the role of gauge invariance as well as the polarization de- 
pendence for anisotropic gaps. In order to compare the 
theoretical predictions with the experiment, we evaluate 
them numerically in a quantitative manner (including ab- 
solute scattering efficiencies!) and compare them to the 
experimental findings. 

Several calculations of the electronic Raman scattering 
efficiency of HTSC have already been published. Some 
of them use highly simplified 2D band structures and 
a decomposition of the Raman vertex 7k in Fermi sur- 
face (FS) harmonicsO or Brillouin zone (BZ) harmon- 
ics, as well as ES integrations instead of the required 
BZ integrationsBEJtj The results of these calculations 
depend very strongly on the number of expansion co- 
efficients piised for 7k and their relative values. Another 
approacht-3 involves the use of band structures calculated, 
in the framework of the local densi ty approximation!!^ 
(LDA) using the LMTO methodML3 Within the approx- 
imations of the LDA, this Raman vertex is exact, i.e. the 
only errors made in such a calculation arise from limita- 
tions of the LDA method itself and from the discretiza- 
tion of the Brillouin zone or Fermi surface. Some of these 
calculations, however, suffer from the fact that only the 
imaginary part of the Tsuneto functioriiia has been used, 
and that only 2D integrations were performed. Q 

The present approacliLj is based on the full 3D LDA- 
LMTO band structure. It uses a BZ integration, screen- 
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ing effects are included, and both the real and imagi- 
nary part of the Tsuneto function are used as required 
by the theory. Electronic Raman spectra are calculated 
for YBa 2 Cu 3 7 (Y-123) and YBa 2 Cu 4 8 (Y-124). The 
orthorhombicity of the cuprates is also taken into account 
in the Raman vertex since we use as starting point the 
band structure of the orthorhombic materials. 

The cuprates under consideration are not only of in- 
terest because of their superconducting, but also of 
their strange normal-conducting properties. Usual met- 
als should show peaks in their Raman spectra at their 
plasma frequencies corresponding to Raman shifts of a 
few eV. The optimally doped cuprates, in contrast, show 
a very broad electronic background (from to about 1 eV 
Raman shift), which is almost independent of temper- 
ature and frequency. The spectra of the underdoped 
HTSC, such as Y-124, show some temperature depen- 
dence at low frequencies {Tilo <C kT). It is possible to ex- 
plain these peculiarities, together with other properties, 
by assuming a certain form of the quasiparticle lifeidxne, 
as was done in the Marginal Fermi Liquid theoryEJEH 

For the superconducting state, various forms for the 
gap function have been proposed. That which has re- 
ceived most experimental support has d x 2_ y 2 symmetry, 
i.e., B\g symmetry in tetragonal HTSC. The power of 
Raman scattering to confirm such gap function has been 
questioned, because, among other difficulties to be dis- 
cussed below, it only probes the absolute value of the gap 
function, i.e. it cannot distinguish between a d x 2_ y 2-\ike 
gap function (for instance cos 20), and a |cos 2<f>\ gap func- 
tion, which corresponds to anisotropic s (A\ g ) symmetry. 
However, it was pointed out that addition of impurities 
can be used to effect the distinction.E3 

This paper is organized as follows: in Sec. E3, we review 
the main properties of the band structures of the investi- 
gated cuprates, as obtained by LDA-LMTO calculations. 
Sec. IH discusses the theory of electronic Raman scatter- 
ing in systems with an anisotropic band structure. We 
first introduce the basic concept of Raman vertex, and 
then present an expression relating the scattering effi- 
ciency to the Raman susceptibility. In Sec. IV and |v| 
we derive expressions for the Raman susceptibility in the 
superconducting and normal conducting phases, and dis- 
cuss some effects not directly contained in the presented 
form of the theory. Section VI is concerned with aspects 
of the experiments which have to be taken care of, espe- 
cially with re gard to the comparison with the theory. Fi- 
nally, in Sec. VII the results of our numerical calculation 



arc presented and compared to the experimental results. 
The difference between calculations involving only FS av- 
eraging, and those in which such averaging is performed 
over the whole BZ, are discussed. 



II. LDA BAND STRUCTURE 

The basis of our calculation is the LDA-LMTO band 
structure of the HTSC under consideration.!!!! For the 
sake of further discussion, we shall describe briefly such 
band structure. _ 

The Fermi surface of YBa 2 Cu 3 7 (Y-123)B consists 
of four sheets, an even and an odd pdcr-like plane band, 
a pda-like chain sheet and a very small pe?7r-like chain 
sheet. The latter is predicted by the full-potential LMTO 
calculations as well as LAPW calculations.!^ We use 
the atomic spheres approximation (ASA) to the LMTO, 
which does not reproduce this rather small feature. In the 
case of Y-123 the three pda-like conduction bands extend 
from — 1 eV to 2 eV relative to the Fermi energy. They 
are embedded in a broad valence band, which ranges 
from — 7cV to 2eV and consists of 36 bands (mainly 
Gu-d and O-p orbitals). Below — 7eV, there is a gap of 
4eV. Above the conduction band, there is another gap 
of 0.5 eV, above which are the lowest fully unoccupied 
bands which consist mainly of d orbitals of Y and Ba. 

The band structure of YBa 2 Cu40g (Y-124) shows sim- 
ilar features. There is an additional pda-\ike chain band, 
while the pci7r-like chain bands are predicted by both, 
full-potential LMTO and LAPW to contain no holes, i.e. 
to be completely filled. 

An interesting feature of the band structure_.of both 
Y-123 and Y-124 is an extended saddle pointO on the 
fo^-axis near the X point. This extended saddle point 
corresponds to a van Hove singularity at approximately 
25meV (Y-123) and 110 me V (Y-124), respectively, be- 
low the Fermi level. As will be shown, the compara- 
tively large density of states in this energy region and 
the warped nature of the corresponding bands has an 
influence on the calculated electronic Raman spectrum. 

The band structure which wc used in our numerical cal- 
culations was evaluated for Y-123 on a mesh of 48 x 48 x 12 
points in the first BZ, involving 4373 irreducible points. 
The band structure of Y-124 is less sensitive to the res- 
olution of the grid (because the extended saddle point 
lies deeper with respect to the Fermi surface). It was 
thus sufficient to use a 24 x 24 x 12 mesh with 1099 ir- 
reducible points. The calculations of the self-consistent 
potential have been performed in the ASA. Therefore, 
the pdir chain band around the S point, which should be 
partly filled, is completely filled. Because of the small 
number of states involved, we do not think this should 
affect significantly our results. 

As stressed above, our calculations are based on a band 
structure obtained within the LDA. We are aware of the 
fact, that the mean free path for transport in the direc- 
tion of the c-axis is smaller than the size of the unit cell, 
i.e. that a description by means of a band structure e n w 
may be questionable (Ioffe-Regel- limit). Nevertheless a 
nontrivial band structure in the c-direction may simulate 
some of the c-direction confinement effects and represent, 
after integration along k z , a reasonable 2D band struc- 



2 



ture. 



III. GENERAL THEORY 

Two approaches have been used to derive the cross 
section (called scattering efficiency when referred to unit 
path length in the solid) of electronic Raman scattering 
in superconductors with anisotrpoic Fermi surfaces. The 
first uses Greenls Junctions ,Bil3EJ and the second the ki- 
netic equation.c3lj Both start with the simplification of 
the Hamiltonian, using kp theory, which relates the. Ra- 
man vertex 7k to the inverse effective mass tensor We 
first briefly review this procedure and, subsequently, the 
derivation of the expression for the scattering efficiency 
using the diagrammatic approach. 

A. The Raman vertex 

To derive an expression for the efficiency for electronic 
Raman scattering, one has to replace the momentum p 
in the Hamiltonian by p — (e/c)A. This yields two dis- 
tinct perturbation terms: Haa = (jo/2)A 2 , quadratic 
in the vector potential, and Ha = — (e/mc)A • p, linear 
in A (we use the transverse gauge; r$ = e 2 /(moc 2 ) de- 
notes the classical electron radius). The relevant states 
in the theory are composed of the state of electrons in 
the sample plus the state of the photon field. The initial 
state of the photon field has laser photons with wave 
vector kj, and polarization and n$ — scattered pho- 
tons with wave vector kg and polarization eg. The final 
state has one laser photon less but one scattered photon 
(Stokes scattering). The vector field can thus be written 
as a superposition of an incoming and a scattered plane 
wave, A = A^ + Ag with 

A s - A+e s e^ r , A L = A^e 1 ^ , (1) 

where Ag contains the creation operator for the scat- 
tered photon and A^ the annihilation operator for the 
laser photons (note that these are not hermitian). 

Since Raman scattering is a second order process in A, 
the term Haa has to be treated in first order perturba- 
tion theory. It is therefore nonresonant and includes only 
intraband scattering. The matrix elements are given by 

M^^^^^iro^k + qlA 2 !^) 

= r Q (A+A L )e s e L S nznf , (2) 

whereas (AgA^ denotes a matrix element involving the 
initial and final state of the photon field and q = k^ — kg 
is the momentum transfer from the photon field to the 
sample. For other values of q, the matrix element (|^) 
vanishes. 

The second term, Ha, produces resonances via second 
order perturbation theory. It has the form 



(q, k) = r < A+A L ) ^U;n m (q, k) (3) 
with the expression 

r( 2 )„ s;nm (q,k)= (4) 
(n f k + q | e* s p | n m k + k L ) (n m k + k L \ e L p | r^k) 

e«ik — En m k+kx, + iO L + iO 

| (nfk + q | e L p | n m k - kg) (n m k - kg | e* s p | njk) 

Eriik — £n m k-k s — Wg + «0 

Here, wl and wg are the frequency of the incoming and 
scattered light, respectively. Note that the states in the 
sum above are states of the sample only. We have used 
Bloch states with band and crystal momentum indices. 
The wavevectors of light, ki and kg, can usually be ne- 
glected in the matrix elements of Eq. (||) because vf <C c. 
For the same reason, e nm .k+k L ~ £ n m ,k- Therefore, we in- 
troduce the symbol Tn f m;n m (k) to denote expression (^) 
with the light wavevectors set equal to zero. 

If we now add the contributions of both terms in 
Eqs. (||) and (|J) and introduce second quantization, wc 
are left with the effective Hamiltonian 

H R = r (A + s A L )~ P(l (5) 

as perturbation Hamiltonian leading to Raman scatter- 
ing. The effective density operator /5 q can be expressed 
in the form 

Pq= 7n/n i (k)c+ ik+q C„ !; ,k ■ (6) 

n i ,ni ,k 

using fcrmionic creation and annihilation operators for 
Bloch electrons as well as the nondiagonal Raman vertex 

lnfHi (k) = e* s e L 5 nfrH + J2 rjg,,^ (k) . (7) 

If we are interested mainly in the low frequency region, 
say Raman shifts below 50 meV, no real interband tran- 
sitions of significant weight are possible. This can eas- 
ily be seen from the band structure (Fig. 2 of Ref. |i"6| ). 
Therefore, we introduce the (intraband) Raman vertex 
7n(k) = 7rm(k). 

We proceed by discussing a very important simplifi- 
cation of ([|) (with m = nf = n and q = 0), the ef- 
fective mass approximation. Four different cases will be 
discussed. First, the virtual intraband transition with 
n m = n i- In this case, up to first order in vp/c, we 
have (n m k | p | n;k) = (njk | p | n m k) (remember that 
n,i — rif) and e ni k — £ n m k = 0. Then, it can be seen that 
the contributions of virtual miraband transitions rela- 
tive to the contribution of virtual interband transitions 
to intermediate states are of the order of the Raman shift 
over the laser frequency, i.e. lo/loj, -C 1, and can therefore 
be neglected. The second case are the virtual mierband 
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transitions involving bands which are much farther away 
from the FS than the light frequency. Then, because of 
\(-m ~ e n m \ 3> k>L, the light frequencies Wi as well as wg 
can be neglected in (|4[). The third case also involves vir- 
tual interband transitions, but for bands at about the 
laser frequency above the Fermi surface. Here, the scat- 
tering is resonant, and the spectra are expected to depend 
strongly on the laser wavelength. One can try to avoid 
this situation by using different laser lines. So we assume 
that in the third case u>l and us also can be neglected. 
Finally, the forth case consists of virtual interband tran- 
sitions to neighboring bands with Ae «C wl. In this case, 
neglecting ujl and u>s is more difficult to justify. We do it 
nevertheless and reach the approximate conclusion that 
we can neglect the light frequencies in Eq. (^) and can 
restrict the sum in (g) to all n m ^ rii. Then, Eq. (|t]) 
becomes completely equivalent to the expression for the 
inverse effective mass from k • p theory and we can write 

T-W^L^^^ (8) 

i.e. the Raman vertex is equal to the inverse effective mass 
contracted with the polarization vectors of the laser light 
and the scattered light, respectively. 

Therefore, using the term Hr with the intraband Ra- 
man vertex 7„k = 7„ n (k) in Eq. (jfy as perturbation to 
the Hamiltonian for A = and treating this in first order 
perturbation theory is, under the mentioned restrictions, 
equivalent-|to taking into consideration both terms H A 
and H AA B 

According to the LMTO calculations, for Y-123 and 
Y-124 there are bands above a band gap between approx- 
imately 2eV and 2.5 eV above the Fermi energy. These 
bands can present a problem with respect to the former 
discussion, because they are almost resonant for typi- 
cal laser wavelengths like 514.5 nm. The same is true 
for the conduction bands, which extend until 2 eV above 
the Fermi surface. Note that due to the strong on-site 
interaction at the Cu-d orbitals, correlation effects are 
expected to be important in the electronic structure. It 
is possible that at energies of the order of 1 eV or more 
above the Fermi surface the picture of the Hubbard bands 
is a better description of the band structure and may ex- 
plain the weak dependence of the Raman spectra on the 
laser frequency observed for laser frequencies in the vis- 
ible range. The band structure shows many bands at 
about the laser frequency below the Fermi energy. These 
should yield resonant contributions to the Raman effi- 
ciency. 

Because the Raman vertex 7k is, in the given approxi- 
mation, the second derivative of the energy with respect 
to k, the Aig component for tetragonal crystals vanishes 
in this version of the theory (^2g is the symmetry of an 
antisymmetric tensor). If one considers once more the 
effects of a nearby resonance, it can be easily seen that 
the Raman tensor does not have to be symmetric. This 



stresses again the questionability of the effective mass 
approach if the scattering is resonant. 

B. The scattering efficiency 

Using the effective mass approach, we arrived at the 
effective Hamiltonian @ with the effective mass deter- 
mining the Raman vertex. This effective Hamiltonian 
is linear in - As- The derivation of the scattering 
efficiency using linear response theory is now a straight- 
forward task. 

The first step is finding a relation between the Raman 
efficiency and a dynamical structure factor of the sample. 
Then, in a next step, the fluctuation-dissipation theorem 
is used to connect the dynamical structure factor to the 
imaginary part of a susceptibility, in our case the Raman 
susceptibility. 

To establish the relation to the dynamical structure 
factor, we add the time evolution factor e~ luJt to the ef- 
fective Hamiltonian (||) and use the golden rule to find 
the transition rate from a state i to a state / of the sam- 
ple. Then, we sum over all final states / of the sample 
and do a thermal averaging over the initial states i. The 
transition rate from a state with til = «k L e L laser pho- 
tons and no scattered photon to a state with ul — 1 laser 
photons and n$ = «k s e s = 1 scattered photon at a tem- 
perature T is given by the expression 

r T (k L ,e L ;k s ,e s ) = ^r 2 • \{A+A L )\ 2 ■ S T (q,u;) (9) 

(the superscript T denotes temperature dependency) 
whereas 

S T (q, W ) = 

E ~z~\(f I p« I *>l s ( E t ~ E i + M ( 10 ) 

i,f 

is a generalized dynamical structure factor (of the sam- 
ple!). The partition function is denoted by Z, and (3 
is the inverse temperature. Now, we sum over all final 
states in a certain region dfl cLus of /c-space around kg 
and normalize to the incoming flux Ticul- This yields the 
expression 

-^-(q,^) = ^r 2 ^(q,c) (11) 

for the differential cross section d 2 a/d£ldu> for a given 
Raman shift lu and a given momentum transfer q. This 
differential cross section is proportional to the scattering 
volume. When performing the calculation for a scatter- 
ing volume equal to unity, a becomes the commonly used 
Raman scattering efficiency. 

Finally, one can define a linear response function, the 
Raman susceptibility 

XRamanM) = ~ ^{Z^e^ [p q (i), p_ q (0)]} (12) 
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and its Fourier-transformed XRaman(q, <*>). To relate the 
imaginary part of this quantity to the structure func- 
tion S T (q, to), we use the fluctuation-dissipation theo- 
rem. The result is 

S T (q,u;) =--(l + fi w )IniXRaman(q,w) (13) 

with the Bose factor n u . 

Equations ( pd| ) and (O) relate the Raman efficiency 
directly to the imaginary part of the Raman susceptibil- 
ity. This evaluation of the Raman susceptibility shall be 
given separately, (i) in Sec. |y| for the superconducting 
phase and Raman shifts of the order of the gap, and (ii) 
in Sec. ^ for large Raman shifts in the superconducting 
phase and for the normal phase. 



IV. THEORY: SUPERCONDUCTING PHASE 

As pointed out in Ref. [lOL the Raman susceptibility 
due to pair-breaking and including screening is given by 
a polarization-like bubble made of a renormalized Ra- 
man vertex Ak, a Raman vertex 7k, and in between 
two Green's function lines for Bogoliubov quasiparticles 
(Fig. |a). The vertex renormalization includes correc- 
tions for Cooper-pair-producing attractive interaction as 
well as the repulsive Coulomb interaction, the Dyson 
equation for the vertex Ak in the limit q — ► is given 
by Fig. 13 in Ref. |(| 

To show more clearly the effect of screening, we write 
the equation for the Raman susceptibility as given in 
Fig. [ljb and 0c. Figure 0b (with a — 7k and b = 7k) 
shows the unscreened susceptibility X77 given by a bare 
polarization bubble with two Raman vertices 7k and the 
contraction of a BCS-like ladder sum with two Raman 
vertices. Therefore, includes the attractive Cooper- 
pair-producing interaction. We include Coulomb screen- 
ing by virtue of a RPA-like sum given in Fig. 0c. The 
effect of screening on the electronic Raman scattering 
can now easily be seen.EI If we denote by Xab a bubble, 
renormalized by pairing interaction, with vertices a and 
b at the ends as in Fig. 0b, the RPA-chain can be easily 
summed up (see Fig. 0c) yielding 

XRaman(q 0, Ul) = % 77 M 1 — — , (14) 

where terms of order q 2 have been dropped. In Eq. ( |T4] ) 
we have used the fact that Vq/(1 — Xn^q) equals 

-Xul 1 - 1 /^. and the factor (1-1/e) is l + 0(q/q T F) 2 - 
Without taking into account Coulomb interaction, 
the Green's functions have a well-known massless pole 
(Goldstone mode) which is a consequence of the break-, 
ing of gauge symmetry in the superconducting phase.c^l 
Coulomb interaction makes this pole acquire a finite mass 
(which can be shown to correspond to the plasma fre- 
quency), so if we correctly include Coulomb screening 



(not done in Ref. £7|) we no longer have a Goldstone 
mode, but a massive Anderson-Bogoliubov mode. This 
mode has the energy h~tu p (u) p is the plasma frequency) at 
the r point and is therefore negligible for the low energy 
behavior of the Raman spectra. 

The susceptibilities \ab in Fig. 0b are like a ladder sum 
contracted with vertices ak and and can be written as 
a sum 

Xafc(q=0, uj) = 2^ a k 6 k A k (w) (15) 

k 

which involves the Tsuncto functionEl Ak(u>). For small 
values of q (compared to the inverse coherence length £ 
and the Fermi wave vector hp), the attractive interaction 
does not have to be taken into account in the summation 
of the ladder, and the Tsuneto function is given simply 
by a unmodified bubble and can be evaluated easily to 
be 

A 2 / i?k \ 
AkM = -g|tanh \^—J x 

X {^2E k + u; + iO + 2Ek-u-io) ' ^ 

Equation ([l6|) involves the gap function Ak (which de- 
pends on the temperature) and the quasiparticle disper- 
sion relation El = £_l + A 2 , with = (e k - e F ) 2 . The 
constants h and ks have been set equal to 1. As al- 
ready mentioned, vertex corrections due to the pairing 
interaction are neglected. This approximation is valid 
for q <C £ _1 ,fci? (Ref. |l^) and uj <C uj p , because the 
Anderson-Bogoliubov pole at the plasma frequency need 
no longer be included. 

A first and very important fact in the expressions above 
is that they contain only the absolute square of the gap 
function, i.e. Raman scattering is not phase sensitive, 
and consequently cannot distinguish between a strongly 
anisotropic s gap |d x 2_j,2| and a d x 2_ y 2 gap. 

In the preceding calculation of the unscreened correla- 
tion functions Xab, we have neglected impurity scattering 
as well as scattering between quasiparticles (collisionless 
regime). In isotropic s-wave superconductors at T = 
and for Raman shifts uj <C 2A, it is perfectly reasonable 
to neglect impurity scattering, because in this regime pair 
breaking is not possible.t3 Also, the scattering between 
quasiparticles can be neglected because their density is 
very small for small temperatures T <C T c . For d-wave 
superconductors this is no longer true. The effect of im- 
purities will be discussed in the next subsection, whereas 
a discussion about scattering between quasiparticles can 
be found in Sec. |v|. 

The second term of (|l4|), representing screening, van- 
ishes if the average of 7k • Ak does. The Tsuneto func- 
tion is fully symmetric, i.e. has A\ g [D^h group) or A g 
(D2h) symmetry regardless of gap symmetry. As a con- 
sequence, the screening term vanishes unless the Raman 
vertex has the same symmetry as the crystal. In the 
tetragonal case. A^-likc vertices are screened, but B\ g - 
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and B2g-like are not. This is different for orthorhom- 
bic HTSC of the YBCO-type. In this case the Tsuneto 
function has A\ symmetry, and the same is true for the 
d a .2_j / 2-like component of the mass [Big of D^h group, A g 
of D2h)- Consequently, in these orthorhombic crystals, 
the Big component is also screened. This discussion is 
also applicable to BISCO, but with interchanged roles of 
Big and B2 g modes because of the different orientation 
of the crystallographic unit cell with respect to the Cu-0 
bonds. 

In tetragonal systems, the Bi g component of the Ra- 
man vertex has nodes at the same position as the gap 
function. This has severe consequences for the low-energy 
part of the spectra.!! 3 ] In two dimensions, the existence of 
the nodes of the gap function in the case of a d x 2_ y 2 gap 
results in a linear density of states at low energies. If 
the vertex has a finite value in this region, the imaginary 
part of the Raman susceptibility is also linear in the fre- 
quency. If the vertex has a node, however, its magnitude 
squared becomes quadratic with respect to the gap on 
the Fermi surface. This causes two additional powers of 
the frequency to appear, the Bi g component of the scat- 
tering efficiency is cubic at low frequencies Two effects 
can alter this behavior: an orthorhombic distortion and 
impurities. 

In our calculations, we focus on a d x 2_ y 2-\ike gap func- 
tion which is only a function of the direction in fc-space, 
but not of the magnitude of k, since the values of the gap 
functions sufficiently far from the Fermi surface do not 
affect the results. We are using the same gap function 
for all bands involved. 



A. Effect of impurities 

In contrast to scattering at non-magnetic impurities in 
conventional (isotropic) superconductors, the influence of 
impurity scattering plays an important role for supercon- 
ductors with anisotropic gaps and its effect on the Ra- 
man spectrum is most pronounced for superconductors 
which exhibit regions in fc-space wher e the gap almost 
or completely vanishes. It was showntBn3 that in the 
case of ci-wave pairing, impurity scattering can be de- 
scribed by extending the nodal points on the 2D FS to 
small finite regions with vanishing gap. This causes a 
nonvanishing density of states at the Fermi energy. For 
anisotropic s-wave pairing the gap anisotropy becomes 
smeared out leading to an increase of the minimum gap 
value A m i n . In the case of a \d x 2_ y 2 1 gap, this minimum 
gap increases monotonically with the impurity concen- 
tration rii for small values of m . 

The renormalization of the gap function by the pres- 
ence of impurities causes an additional contribution, 
which is linear in the Raman shift u for small Raman 
shifts u>, in the Raman spectra.E3 This has consequences 
for the Big spectrum of a tetragonal crystal, which, ac- 
cording to the theory, has a cubic ^-dependence, because 



a linear frequency dependence is added. As will be dis- 
cussed in the next subsection, the orthorhombicity of the 
YBCO compounds also causes a linear addition to the 
cubic behavior of the Bi g channel spectrum. 

In the case of a \d x 2_ y 2 |-like, A g symmetry gap func- 
tion the impurity-induced minimal gap A m i n causes an 
excitation-free region to show up in the electronic Raman 
spectrum below a Raman shift of 2A m i n . 



B. Effect of orthorhombic distortion 



As already mentioned, orthorhombic distortions, i.e. 
deviations from the tetragonal symmetry, have a dif- 
ferent effect on Y-I23 and on Bi-22I2. Consider the 
Big {D^h) component of the inverse mass tensor in a 
tetragonal high-T c superconductor with a dj^^-like gap. 
The Big (D4h) mass has its nodes in directions diagonal 
to the axes of the copper planes; the same is true for 
the gap function. As mentioned above, this results in 
the w 3 -dependence of the Raman efficiency for Bi g (-D^) 
scattering, in contrast to the ^-dependence predicted for 
Ai g and £?2g scattering. Let us now consider the or- 
thorhombic distortion present in Y-123. The zeros of the 
Big {Dnh) mass shift because there are no longer mirror 
planes through the (110) axes. For this reason, the low- 
energy part of the spectrum acquires a linear component 
in addition to the w 3 component of the D^h case. 

In Bi-2212 the situation is different because the or- 
thorhombic crystallographic cell is rotated by 45° with 
respect to the a- and 6-axes: the orthorhombic distor- 
tion preserves the mirror planes [a ± b, c] . Consequently, 
the Bi g zeros stay at the same position, the low-energy 
efficiency acquires no linear component. 



C. Effect of multilayers 

In systems with one layer of CU-O2 planes per unit 
cell there is only one sheet of Fermi surface and the 
mass fluctuations are essentially intraband mass fluctua- 
tions, which are very sensitive to the scattering polariza- 
tions. The scattering related to the average mass is fully 
screened. The simplest A\ g (-D47J scattering is related 
to a Raman vertex of the form cos 40 symmetry while 
Big (Dih) scattering is obtained for a cos 20 vertex. In 
multilayer systems, interband fluctuations between the 
various sheets FS are also important. The lowest compo- 
nent of such fluctuations corresponds to different average 
masses in each FS sheet. Such fluctuations do not depend 
on the scattering polarizations and lead to unscreened 
scattering of A g symmetry. 



G 



D. Effect of sign change of 7^ on the Fermi surface 

The behavior of the Raman vertex near the Fermi sur- 
face, especially its sign, is crucial for the scattering effi- 
ciency and, in particular, for the effect of screening. Anti- 
screening, i.e. an enhancement of the scattering efficiency 
by screening, can occur if the Raman vertex changes sign 
on the Fermi surface. This can be seen by considering the 
screening part 



Im XSc 



Im 



X71 



(17) 



of the Raman susceptibility. A positive value of Im xscr 
enhances the efficiency, i.e. corresponds to antiscreening. 

To show how antiscreening arises, we first write the 
screening term Im Xscr in terms of the real and imaginary 
parts A' = Re A and A" = Im A of the Tsuneto function 
and the Raman vertex 7 as 



Im XSci = 



( 7 A') 2 (A")-(7A") 2 (A") 



2( 7 A')(7A")(A') 



(AT + (A'T 



(18) 



The imaginary part of the Tsuneto function A" is a posi- 
tive (5-function. Consequently, the quantity (A") is a pos- 
itive function of the Raman shift uj. If 7k changes sign 
in a region around the Fermi surface, it is possible that 
(7A") changes sign as a function of ui, i.e. has a zero. At 
the position of this zero, the second and the third term in 



the numerator of (18) vanish. The first term, 



hat 

is positive and can become dominant in Eq. ( |18| ) . In this 
case antiscreening results. In the Appendix A will be 
shown that antiscreening is particularly sensitive to the 
sign of the Raman vertex on parts of the Fermi surface 
around the directions of the nodes of the gap function 



V. THEORY: NORMAL PHASE 

In the normal phase, the exact mechanism which 
produces a finite Raman intensity almost constant 
over a broad frequency and temperature range, is 
not known. Therefore, we assume some scatter- 
ing mechanism, which implies a finite lifetime of the 
quasiparticles. Candidates for this scattering are the 
quasiparticle-ouasiparticle scattering in Marginal Fermi 
Liquid theory^ (MFL), impucity scattering^ or scatter- 
ing due to spin fluctuations EH A self energy with non- 
vanishing imaginary part yields a susceptibility of the 
form 



Xab(q=0,w) = y^ak&k^kM 



(19) 



with the relaxation kernel (the function /' is the deriva- 
tive of the Fermi function with respect to the energy) 



= -/'(&) 



ui + ?T k 



and its imaginary part 

Im^(w) = -/'(&)- 



> 2 + rt 



(20) 



(21) 



This can easily be seen by evaluating a bubble with two 
Greens function lines for quasiparticles with an imagi- 
nary part r k of the self energy. 

Note that in the superconducting phase for Raman 
shifts larger than ~ A, the relaxation effects described 
by ( |2~l| ) are also of importance. The relevant relaxation 
kernel in this case is 



^kM = ~f'(E k ) 



iTi, 



iTi, 



(22) 



where £ k = (e k - e F ) 2 . 

To describe the constant background in the Raman 
spectra in the normal phase, one has to adooUthe quasi- 
particle scattering rate of the MFL theoryBEll 



r k (w) ~ max(aT, (3ui) . 



(23) 



In order to evaluate the real part of j/ k using causality 
arguments, and to prevent divergences, we introduce a 
high-frequency cutoff ctto_Note that the nearly antiferro- 
magnetic Fepni liquicOo (NAFL) and also the nested 
Fermi liquid^ (NFL) yield a very similar quasiparticle 
scattering rate. The former can also provide a mecha- 
nism, which accounts for d x 2_ y 2 pairing. Similar-Results 
are obtained with Luttinger liquid based results. E^l 

Equation (|2^) yields a scattering continuum which is 
constant for frequencies smaller than min(aT//3, T) and 
for frequencies larger than the temperature T, but with 
different intensities. In the first case, r k is proportional 
to the temperature, i.e. Imx ~ w/T. Multiplying by the 
Bose factor 1 + n w ~ T/uj a constant is found. In the 
second case, r k ~ uj, and, consequently, Imx = const. 
The Bose factor is also constant and one is left with a 
constant Raman intensity. Note that in the first case, 
Imx cancels the ui- and Independence of the Bose fac- 
tor. It has been shown,L3E30 that YBa2Cu40s does 
not exhibit this behavior. This has been attributed to 
the breakdown of MFL theory for not optimally doped 
cuprates.EJ Actually, in this case the spectra are nearly 
temperature independent after dividing them by the Bose 
factor. We shall address this question once more at the 
end of this section. 

To discuss quasiparticle-quasiparticle (qp-qp) scatter- 
ing, and its influence on electronic Raman scattering, we 
start with the case of a d x 2_ y 2 gap. Suppose the nodes 
of this gap have a width Sq in fc-space on the Fermi sur- 
face due to impurity scattering. We use the model of 
Eq. ( p2[ ) with a quasiparticle scattering rate T k indepen- 
dent of k and discuss first the case T — 0. Then it can 
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be seen that the contribution of qp-qp scattering to the 
imaginary part of the Raman susceptibility (|l^) for low 
frequencies lu <g; A max is proportional to the Drude-like 
factor luT / (w 2 +r 2 ) (which is, for small lu and low temper- 
atures T < ui, linear in lu if T — const (semiconductors) 
or r ~ max(o; 2 ,r 2 ) (FL), but constant as a function of 
lu if r ~ max(cj,T) (MFL). In the tetragonal case, it is 
also proportional to the density of states at the Fermi 
surface and in the case of Ai g and Big polarizations to 
the width Sq, and in the case of Bi g to the third power <5q 
of the width So. The discussion for BISCO is analogous 
with the exception that B\ g and Big exchange their role. 

Finite, but small temperatures T <C A max have the 
effect of enlarging the widths Sq linearly in temperature, 
i.e. the temperature dependence of the contribution from 
qp-qp scattering is proportional to const + T. Note that 
for T ;> 0, the Bose factor changes the linear-in-w depen- 
dence to a constant. 



For the anisotropic s gap of the form \d x 



which 



acquires a finite minimum gap A m i n due to the presence 
of impurities,E3 the situation is different. The frequency 
dependence is also given by the factor luT /(lu 2 + T 2 ) in 
addition to the Bose factor. But the temperature depen- 
dence is different. For temperatures T <C A m j n smaller 
than the minimal gap, the density of quasiparticles is 
proportional to exp(— A m j n /fcT), i.e. the contribution of 
qp-qp scattering to the Raman efficiency is exponentially 
small. At kT « A m ; n , this exponential dependence on T 
crosses over to a power law. 

The background electronic Raman spectrum in the 
normal phase is almost independent of temperature for 
nearly optimally doped high-T c compounds only. In the 
overdoped and underdoped case, the materials seem to 
show Fermi liquid-like behavior concerning the quasipar- 
ticle scattering rate Tk (for small w).e1tE3 The temper- 
ature dependence of the scattering rate Tk, as defined 
in (^2|), has been measuredEH for optimally doped and 
overdoped Bi-2212, and, especially in the case of the 
Big (D^h) mode, the optimally doped sample shows V = 
aT, whereas for the overdoped sample T = a'T 2 + Tq. 
Therefore, the overdoped sample shows properties of a 
normal Fermi liquid which are predicted by theory to 
have r ~ max(w 2 ,T 2 ). The B\ g (D^h) mode result for 
the optimally doped sample yields the puzzling quasipar- 
ticle scattering rate T = const. 



VI. EXPERIMENTAL SPECTRA 

The experimental determination of absolute Raman 
scattering intensities is plagued by a number of difficul- 
ties (a reason why usually "relative units" are found in 
the literature). The first is related to the presence of 
elastically scattered light in the spectra, in particular 
when non-ideal sample surfaces are involved. Depend- 
ing on the quality of the spectrometer this leads to con- 
tributions extending typically, for the parameters of the 



present work, up to 50 cm -1 from the center of the laser 
line. These contributions can be filtered out using a pre- 
monochromator or notch filters but, in any case, Raman 
scattering measurements below 50 cm -1 remain difficult. 
The measurements discussed here have been performed 
by comparison with the known efficiency of silicon after 
correcting for differences in the scattering volumes. The 
procedure leads to errors of about 50%. 

We use for comparison with, the calculation the exper- 
imental data qf-Krantz et al.a in the case of Y-I23, and 
Donovan et alSB in the case of Y-I24. Our Figs. | and | 
are taken from these publications. In the case of Fig. 4 
we have corrected a scale error in the abscissa found in 
Ref. ||. In the case of Fig. || we have calculated the Ai g 
component from the experimental results for the (x'x') 
and (xy) polarizations. 

The classification of the measured spectra according to 
irreducible representations of the symmetry group of the 
crystal is performed with the use of the Raman tensor R 
which is related to the Raman efficiency through the ex- 

~ 2 

pression / ~ |e£i?es| , bilinear in the Raman tensor. In 
the calculations, the Raman tensor does not appear ex- 
plicitly, the inverse effective mass d 2 E / '(dkidkj) playing 
its role. It is important to note that the Raman efficiency 
as given by the theory (Eqs. (|lf|), (|I|), and @) is bi- 
linear in the inverse effective mass of the Raman vertex 
(including the screening part!), i.e. contains the same in- 
terferences as the approach involving the Raman tensor. 
Note that the Tsuneto function A is fully symmetric. In 
the normal phase, the scattering kernel v has been as- 
sumed to be the same for all scattering channels. 

In most of the measurements of the Raman efficiency 
in orthorhombic high-T c superconductors, an A\ g com- 
ponent has been given. Strictly, this irreducible repre- 
sentation does not exist in Dih but only in D±h. In 
orthorhombic crystals, the Raman tensor contains two 
A g components which correspond to the A\ g and B\ g 
components of the tetragonal D±h case, and which are 
not distinguishable in Dih because they transform in the 
same way. Nevertheless, quantities can be constructed in 
the orthorhombic case which correspond to the tetrago- 
nal Ai g component. 

One of these is 1^ = (I xx + l yy )/2 — I x i y >. Both, I xx 
and Iy V contain A\ g and B\ g (D^h), and also an inter- 
ference term which cancels when I xx and I yy are added. 
The I x 'y' efficiency contains B\ g and Ai g (D^). If we as- 
sume that the antisymmetric component (A^g in D^h) of 
the Raman tensor R vanishes (i.e. I xy — I yx ), I x > y ' corre- 
sponds to tetragonal B\ g and cancels the B\ g contribu- 
tion in I xx and I yy . Provided that the Ai g component of 
the Raman tensor vanishes, 1^ corresponds to the Ia 1b 
of the tetragonal case. Note that the antisymmetric com- 
poment (R xy — R yx )/2 of the Raman tensor vanishes in 



the effective mass vertex theory given in Sec. Ill because 
of 

Ixy — lyx regardless of the symmetry of the crystal, 
and also in the experiment in the case of tetragonal crys- 
tals but not necessarily for orthorhombic crystals. The 
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equality of I xy and I yx in the calculation is an artifact of 
the theory. 

A second possible construction for A\ g is I*- 2 -* = I x > x > — 
I xy . The I x ' x i efficiency contains A\ g and Bi g contribu- 
tions. The interference term of these two contributions 
vanishes in the tetragonal as well as the orthorhombic 
case. Both, Bi g (D^h) and A^ g are contained in I xy . But 
if the Ai g component of the Raman tensor vanishes, 1^ 
also corresponds to the lA n of the tetragonal case. In one 
of the experimental worksa a different method to extract 
the A\ g component was used. Both of the expressions for 
and /( 2 ) contain contributions of the A2 g (D^h) Ra- 
man tensor component. This component may be present 
in the experiment, but not in the theory, a fact, that 
has to be kept in mind when comparing the numerical 
results to the measurements. Note that the Raman ef- 
ficiencies in (xy) and (x'y 1 ) polarization configurations 
also contain contributions from the antisymmetric part 
of the Raman tensor. In view of these uncertainties in 
A\ g we mainly focus in the next section on the directly 
observable components of the Raman tensor. 

We shall conclude this section by taking up again the 
question of the validity of the effective mass approxi- 
mation. In the experiment, this can be checked in two 
ways. First, via the dependence of the spectra on the 
laser frequency which should make it possible to distin- 
guish the contributions to the Raman efficiency resulting 
from resonant and non-resonant transitions, respectively. 
The second way involves the measurement of the A2 g 
component of the mass. If the effective mass approxi- 
mation is valid, the Raman vertex should be symmetric 
{"Ixy = lyx), i.e. the A2 g (D±h) component should van- 
ishes. A non-vanishing A2 g component of the measured 
scattering would cast doubts on the appropriateness of 
the effective mass approximation. 



VII. NUMERICAL RESULTS AND DISCUSSION 

To carry out the numerical BZ and FS. integrations, 
we employed a tetrahedron approach.HEj The conver- 
gence of the integrations was checked by using different 
meshes. In Figs. ^ and |^, the results of full BZ integra- 
tions for Y-123 and Y-124, respectively, are plotted. The 
corresponding spectra obtained through FS integrations 
can be seen in Ref. [51]. The Bose factor has not been in- 
cluded, hence the results apply to zero temperature. In 
both figures, the Raman shift is given in units of the gap 
amplitude Ao . Since the calculated scattering efficiencies 
for BZ integrations, contrary to FS integrations, are not 
only a function of the reduced frequency but depend also 
weakly on the value of Ao, we took for the calculations 
Aq = 220 cm^ 1 . This value of Ao falls in the range of 
Ao's determined by Raman scattering and other meth- 
ods. The delta-function peaks in the Tsuneto function 
have been broadened phenomenologically by introducing 
a finite imaginary part T — 0.3Aq of the frequency vari- 



able LO. 

Figures || and || display spectra for each of the polar- 
ization configurations (yy), {x'x'), (xx), {x'y'), and (xy), 
as well as the symmetry component A\ g (D^h) (defined 
by Ia 1s — I x > x > — I X y), the unscreened intensities, the 
screening part (jl7|), and the total intensities, equal to 
the difference between unscreened and screening parts. 
Note that the (x'y') configuration corresponds to the 
Bi g (D^h) component because of the vanishing of the A2 g 
component in the theory. 

We discuss first the results for Y-123. The A\ g compo- 
nent (in the rest of this section we use tetragonal nota- 
tion unless explicitly stated) is subject to rather strong 
screening, however its unscreened part is comparable to 
that of the B\ g component. The relation between the un- 
screened and the screened (total) spectral weight of the 
A\ g component is about three. Nevertheless, the shapes 
of the unscreened and the screened parts are the same 
and, consequently, there is almost no shift in the peak po- 
sition due to screening (contrary to the results of Ref. ||). 
The peak is located almost exactly at 2Ao. Note that 
there is no antiscreening in the A\ g component. The low- 
energy part of all A\ g spectra (screened and unscreened) 
is linear, as predicted by the theory. 

As already mentioned, the (x'y') component (equal to 
the Bi g component in the non-resonant case) is almost 
four times stronger than its screened Ai g counterpart. 
The screening is very small, its nonvanishing being an 
effect of the distorted tetragonality of the crystal. There 
is, in this case, a very small amount of antiscreening in 
the region below 2Ao. As in the case of the A\ g com- 
ponent, the (x'y 1 ) component peaks at almost exactly 
the 2Ao frequency shift. The low- frequency part has an 
auj + /3L0 3, frequency dependence, the linear part arising 
from the distorted tetragonality, i.e. the fact that the B\ g 
mass does not vanish at exactly the same position on the 
Fermi surface as the gap function does. 

The efficiency of the peak in the (xy) configuration 
(equal to the B2 g component in the non-resonant case) is 
also four or five times smaller than that of the A\ g peak. 
The (xy) peak is located at about 1.3Ao, as expected 
from the fact that in the neighborhood of the region 
where the gap is large, the i?2g mass vanishes. Conse- 
quently, the peak is not as sharp as in the former cases 
and screening vanishes since these spectra correspond to 
a nonsymmetric (B\ g ) representation of the orthorhom- 
bic group (D 2 h)- 

In the A\ g and (x'y 1 ) spectra there should be a small 
peak at about w — 2-\/e 2 H + A 2 iax « 3.9A due to the 
van Hove singularity on the fc^-axis near the X point. 
The corresponding structure, however, is very weak, and 
practically invisible in Fig. 0. This is not unexpected for 
a 3D calculation. These peaks appear strongly when-2D 
calculations are performed through BZ integrations £3 

In general, the efficiencies in Y-124 (Fig. ||) are about a 
factor of three less than those for its Y-123 counterpart. 
Moreover, the screening of the Ai g component of Y-123 







is much stronger than that of Y-124. This may be, at 
least in part, due to the additional chain band: The (yy) 
component of Y-124 is less screened than the (yy) com- 
ponent of Y-123. At low frequencies, we correspondingly 
have antiscreening even in A\ g , a fact which reveals a 
change of sign of the effective mass on the Fermi surface 
(see Sec. IV D). Due to this antiscreening, the peak in 
the Ai g spectrum is shifted from 2Ao towards approx- 
imately I.6A0. In contrast to the situation in Y-123, 
the Y-124 spectra show clearly the influence of the van 
Hove singularity on the spectra, as a small hump (vH) 
located near 2^/e^ H + A^ ax « 7A . In the A lg spectrum 
this hump is almost screened out whereas in the (x'y 1 ) 
spectrum it appears slightly increased by the influence of 
antiscreening. 

To compare these predictions with the experiment let 
us first focus on the peak positions. The experimental 
results for Y-123 (Fig. |J, lower part) clearly show that 
the position of the (yy), (x'x') and (xx) peaks is at about 
300 cm -1 , whereas the (x'y 1 ) peak is located at 600 cm -1 , 
i.e. at twice the frequency of the former. This fact is in 
sharp contrast with the calculated spectra and has been 
at the-ceater of the controversy concerning the topif 
handO'Ej It has been suggested by Devereaux et 
that the B\ g component peaks at 2Ao, and the A\ g com- 
ponent becomes shifted down to almost Ao by the screen- 
ing. This interpretation contradicts our numerical results 
which clearly suggest that the influence of screening on 
the position of the A\ g mode is usually smaller. The 
frequency renormalizations of phonons around T c also 
seem to contradict .-the interpretation in Refs. || and [l^. 
It has been showrcd that lowering the temperature of 
the sample in the superconducting phase causes the Ai g 
435 cm -1 phonon (plane-oxygen, in-phase) to shift up in 
frequency and the Bi g (D^ notation) 340 cm -1 phonon 
(plane-oxygen, out-of-phase) to shift down. This, in turn, 
implies an amplitude of the gap 2Ao between 300 cm -1 
and 360 cm -1 and is consistent with our interpretation of 
the electronic Raman spectra with the A\ g peak at 2Ao. 

Note that the (yy), (x'x') and (xx) spectra do not con- 
tain contributions of the A^ g (D±h) antisymmetric com- 
ponent of the Raman tensor while the (x'y') component 
does. So, the experimental results may suggest that the 
shift of the position of the (x'y') spectrum with respect 
to the peak position of the other spectra is due to res- 
onance effects. The (xy) spectrum is also influenced by 
the A-i g component. It is difficult to determine its peak 
position from Fig. |^, but it seems to be located at the 
same position as that of the (yy), (x'x') and (xx) con- 
figurations. The calculation predicts it to be located at 
about 1.3Ao, the shift to 2Ao can also be attributed to 
the existence of an Ai g component, like in the case of the 
(x'y') configuration. 

To compare the relative intensities of the spectra with 
different polarizations, we refer to Table |J, which lists 
them together with the corresponding absolute intensi- 
ties, both at the peak position. The detailed results of 



our FS integration have already been reported earlier .El 
We begin with Y-123 (upper panel in Table ||) and com- 
pare BZ integration results to the experimental ones. 
With the possible exception of the A\ g component (and 
the (x'x') component, which is very similar to A\ g ), the 
agreement is rather good. The deviation of the A\ g com- 
ponent may be attributed to screening, which is very sen- 
sitive to sign changes and other details of the Raman ver- 
tex near the Fermi surface (such as details of the band 
structure and especially the exact position of the Fermi 
energy). 

The second compound, Y-124 (lower panel in Table ||), 
also shows reasonable agreement between the results of 
the BZ integration and the experiment. However, we 
also have problems with the A\ g component, as we did 
for Y-123. 

The measured absolute intensities agree particularly 
well with the calculations in the case of Y-123. With the 
exception of A\ g , the discrepancy between theory and 
experiment is only a factor of two, which can easily be 
related to the difficulties in measuring absolute scatter- 
ing cross sections. In the case of Y-124, the discrepancy 
is a bit larger, but a factor of four can still be considered 
good. We should also keep in mind that resonances of ojl 
or lus with virtual interband transitions are expected to 
enhance the simple effective mass Raman vertex, a fact 
which could also explain why the measured scattering 
efficiencies are usually larger than the calculated ones. 

We close the discussion of the numerical results with a 
remark about the Fermi surface integration. For Y-124, 
the results of the former correspond rather closely to the 
results from the BZ integration. The situation is different 
for Y-123. Here, the (xx) peak height is almost a factor 
four larger in the FS integration than in the BZ integra- 
tion. This is likely to result from the close proximity of 
the van Hove singularity to the FS in the case of Y-123 
(25meV), as compared to Y-124 (110 meV). 

To verify the predictions related to the effec t of or- 
thorhombic distortions as discussed in Sec. IV B, we per- 
formed a fit of the function aLU+(3u> 3 to the low- frequency 
part of the B\ g data for Y-123 reported in Ref. y and 
Ref.|42]as well as for Bi-2212 (taken from Ref. f|| and to 
the results of our numerical calculations for Y-123. The 
ratios of the cubic vs. the linear part (at u> — 300 cm -1 ) 
of the fit to the low-frequency efficiency are given in Ta- 
ble |§ 

Both measurements for Y-123 agree in their large lin- 
ear part, which should be due mainly to the lack of exact 
tctragonality and the presence of impurities. The results 
of the BZ integration show a smaller linear part, because 
they do not take into account the influence of impurities. 
Finally, the result for Bi-2212 is completely different from 
the former results for Y-123. The linear part almost van- 
ishes, in agreement with the preceding discussion. 
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VIII. CONCLUSIONS 



In spite of the striking ability to predict not only gen- 
eral features of the observed spectra but also their peak 
intensities, our calculations are not able to predict the 
relative positions of the A\ g and B\ g peaks. Accord- 
ing to Figs. H and || the A\ g spectrum should peak only 
slightly below 2Ao while Bin should peak at 2Ao. The 
experimental data of Figs. [| and g, however, indicate 
that the B± g spectra peak nearly at twice the frequency of 
A\ g . Since the observed A\ g peak is considerably sharper 
than that of B\ gi we may want to assign the A\ g peak to 
2Ao. Our calculations show that it is impossible to repro- 
duce both peak frequencies with a simple gap of the form 
Ao cos 2(j) where </> is the direction of the k- vector. A rea- 
sonable fit was obtained in Ref. |^ with a two-dimensional 
FS which did not take into account the chain component 
and assigned d- and s-like gaps to the two bonding and 
antibonding sheets of the FS of the two planes in an ad 
hoc way. Within the present 3-dimensional band struc- 
ture the FS cannot be broken up into bonding and anti- 
bonding plane and chain components since such sheets 
are interconnected at general points of fc-space. It is 
nevertheless clear that there is no reason why the gap 
function should be the same in the various sheets for a 
given k-direction. Thus the remaining discrepancy in the 
peak positions between theory and experiment could be 
due to a more complicated A„k than a simple Ao cos 20 
used here. Another possible source of this discrepancy 
is scattering through additional excitations of a type not 
considered here (e.g. magnetic excitations) contributing 
to and broadening the B lg peak. 

A BCS-like theory, which involves an attractive pairing 
potential as well as the repulsive Coulomb potential and 
uses an anisotropic d x 2_ y 2-\ike gap function in connection 
with the effective mass approximation used in the calcu- 
lation of the absolute Raman scattering efficiencies yields 
result which are in significant agreement with the experi- 
mental spectra. One exception, the peak positions of the 
Ai g and the B\ g components, remains unexplained. The 
theory predicts them to be both located near u> = 2Ao, 
but the experiment shows the peak in B\ g at almost twice 
the frequency of the peak in A\ g . The weak Bi g spec- 
trum agrees in intensity and peak position with calcula- 
tions for a d x 2_j / 2-like gap. The results of other experi- 
ments, invoking the temperature dependence of phonon 
frequencies,!^ suggest that the A\ g peak position corre- 
sponds to the gap amplitude 2Ao. The shifting of the B\ g 
peak towards higher frequencies may have an origin dif- 
ferent from the mass-fluctuation-modified charge-density 
excitations described in the theoretical part of this paper 
but could also be due to a multi-sheeted gap function, 
more complicated than the simple d x 2_ y 2-Y\ke Aocos20 
gap assumed in our calculations. The initial variation of 
the A\ g and B\ g scattering efficiencies vs. oj are linear as 
expected for that gap. The B\ g symmetry becomes A g in 
the presence of the orthorhombic distortion related to the 



chains. Consequently, the scattering efficiency at low fre- 
quency is not proportional to to 3 but should have a small 
linear component which is found both in the calculated 
and the measured spectra. In the corresponding spec- 
trum of Bi-2212, with and orthorhombic distortion along 
(x + y), the Bi g (D^h) excitations also have a nonsym- 
metric B\ g {Din) orthorhombic character. Consequently, 
for small u> no component linear in u> is found in the mea- 
sured spectra. 

We have performed our calculations using either BZ 
or FS integration. In the case of Y-124 the spectra so 
obtained are very similar. For Y-123 quantitative differ- 
ences appear; they are probably related to the presence 
of a van Hove singularity close to the FS. These singular- 
ities appear as weak structures in the calculated spectra, 
as expected for a 3D band structure. 
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APPENDIX A: ANTISCREENING AND THE 
SIGN OF THE RAMAN VERTEX 



In Sec. IV D. we pointed out that the effect of antis- 



creening results from sign changes of the Raman vertex 
7k, i.e. the inverse effective mass, on the Fermi surface. In 
this appendix we present a different proof using a power 
expansion of 7k- 

For very low frequencies lo <C Ao, only the regions 
around the node directions of the gap function (we as- 
sume a d x 2_ y 2-like gap) contribute to the susceptibility. 
We focus on a specific node of the gap function and 
define the point ko as the point of intersection of the 
node line of the gap function and a specific sheet of the 
Fermi surface. Then we introduce an orthogonal coor- 
dinate system {k,k±^} in the a-&-plane in fc-space with 
the origin at ko, rotated in such a way that the k±- 
axis is perpendicular to the node line (i.e. tangent to 
the Fermi surface). We write the Raman vertex as a se- 
ries 7(k) = J2iJi(k±/k c y (where fc c is a cutoff), using 
the assumption that d^f/dk — 0. This approximation 
is justified since contributions to the Raman susceptibil- 
ity mostly arise from a narrow region around the Fermi 
surface). We write Eq. (J18|) as 



Im XSCT = £lm X « 



(Al) 



(i) 

whereas Im Xsci contains i-th powers in k± from the ex- 
pansion of 7k- The first three terms in this sum are 



11 



Im XScr 

T (!) 
Im XScr 

T ( 2 ) 



-7o 2 (A") 

-2 7o7l (fc ± A") 

-2 70 72<fciA") 



(A2) 



"7r 



(\") ({k ± \") 2 - (k±X') 2 ) + 2(\>)(k ± X)(k ± \") 
(A') 2 + (A"} 2 



Now we investigate the behavior of (fc^A") in the low 
frequency limit. We set vf = 1 and Ao = 1 (this repre- 
sents a simple change in scales). Then, E^ = k (E^ is 
constant as function of k± by definition), Ak = k±, and 
therefore AjV^ = k\ /k 2 and A" - (k ± /k) 2 ■ S(2k - ui) 
from Eq. (|l6|). We perform a 2D BZ integration which 
has to be cut off at a value proportional to uj in the k± 
integration and find 



k 2 

(k l j_X") ~ / dk ± dk k\ ■ -AS(2k - ui) 



dk± k 



: -+ 2 ~ r 



(A3) 



as the low frequency behavior. The same is true for the 
real part \kj_X'). Note that the proportionality constant 



in (A3) is positive. This implies that Imxsc r ~ and 



+0(c 4 )) 



(^27072 + a 3 j 2 )u> 3 
(A4) 



with positive constants ao, a\ and a2- The sign of a 3 
depends on the specific case. 

For a tetragonal Fermi surface, we only have to focus 
on the A\ g mode because the screening contributions to 
the other components vanish by symmetry. The A\ g sym- 
metry implies 71 = because of the ad symmetry opera- 
tion (reflection at TS) which transforms 7 ifcj_ — > — 7ifej_. 
Then, the screening term for low frequencies can be writ- 
ten as 



(a 7 2 w + a 2 7o72^ 3 H ) 



(A5) 



If there are nodes in the Raman vertex near the node 
of the gap function, then 7072 < 0, i.e. the screening 
term is negative for very small ui, but eventually crosses 
zero because of the uj 3 contribution. For large uj <; 2Ao, 

the approximation (7A) « ( 7 )(A) yields Xscr = -(7) 2 Xn, 
which is the screening term for an isotropic Fermi surface 
with a scalar Raman vertex (7}, and therefore negative. 
If the Raman vertex does not show nodes near the node 
of the gap function, 7072 is larger than zero, and the 
screening term is not likely to change sign. 

In the orthorhombic case of a weakly distorted tetrago- 
nality, the Bi g (D^h) zero in the Raman vertex may shift 
with respect to the gap node. This can be described by 
a small 70 ^ and 71 7^ 0; 72 can be neglected. The low 
frequency screening term then has the form 



The first an the third term have almost A\ g symmetry 
whereas the second term has almost B\ g [D^h) symmetry. 
It can be shown that the term proportional to ui 2 van- 
ishes if the other nodes of the gap function are taken into 
account. Whether antiscreening exists or not depends on 
the sign of a 3 . 

It can be seen that for small values of 70, the antis- 
creening can start already at very small Raman shifts as 
is the case in the calculations, Fig. ||[ B\ g panel. 



Im xsc 



(a 7 2 w + ai7o7i^ 2 + a 3 ^\uj 3 ■ ■ •) . (A6) 
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TABLE I. The ratio between the linear and cubic parts of the low energy Raman efficiency in B\ g (D^h) configuration of 
several high-T c compounds at a Raman shift of lo — 300 cm -1 . 



HTSC 


cubic:linear in B\ g 


Reference m 


Y-123 


1 


Krantz et aM 


Y-123 


1 


Hackl et alB 


Y-123 


0.35 


BZ integrati(w*i 


Bi-2212 


0.07 


Staufer et allB 



TABLE II. Comparison of the experimental peak scattering efficiencies given in units of 10 8 cm cm 1 sr 1 to the theoretical 
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0.24 
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0.54 


12 


0.30 


A lg 
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0.96 
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2.3 
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0.31 
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1.1 


0.17 


1.0 


0.23 
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FIG. 1. Incorporation of screening effects into the theory of Raman scattering by electronic excitations in HTSC. The grey 
shaded bubbles are sums of ladders contracted with vertices a and b. Wavy lines correspond to the long-range Coulomb 
interaction, dashed lines to the attractive pairing interaction. The equation on the last line corresponds to Eq. (114). 



FIG. 2. Results from the BZ integration for Y-123. Given in the five panels are absolute efficiencies for electronic Raman 
scattering. The upper three curves are labelled using the irreducible representation (D^) of the scattering mass, the lower 
two panels with the polarization geometry. Each of the five panels contains the total absolute Raman efficiency according to 
Eq. (JTTJ) and (jl3|) and its two constituents, the unscreened and the screening part according to Eq. (p"4j). 

FIG. 3. Results from the Brillouin zone integration for Y-124. For details see the caption of Fig. ^| 

FIG. 4. Experimental Raman scattering efficiencies for Y-123 from Ref. ^j. The vertical scales are absolute Raman efficiencies, 
measured at T — 10 K and an exciting laser wavelength of A = 488 nm (Note that a scale error found in Ref. |^ has been 
corrected). The A\ g component extracted according to Ia\„ = (Ixx + I yy )/2 — I x i y i is plotted in the lower panel together with 
the quasitetragonal B\ g and B29 components. 



FIG. 5. Upper panel: experimental absolute Raman efficiencies given for the five specified polarization configurations 
for Y-124 from Ref. |35| These data are taken at T = 10 K with an exciting laser wavelength of A = 514.5 nm. Lower 
panel: smoothed curves and the Ai g spectrum additionally extracted from the former. In both panels, consecutive offsets of 
0, 1, . . . , 4 x 2.5 cm - sr~ cm were used. 
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Figure 3 
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Figure 4 
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